High resolution amplitude and phase gratings in atom optics 
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An atom-field geometry is chosen in which an atomic beam traverses a field interaction zone 
consisting of three fields, one having frequency = c/A propagating in the z direction and the other 
two having frequencies f2 + Si and fi + 82 propagating in the -i direction. For ni<5i + 712^2 = and 
\Si \ T, \S2\T 3> 1, where m and ri2 are positive integers and T is the pulse duration in the atomic 
rest frame, the atom-field interaction results in the creation of atom amplitude and phase gratings 
having period A/[2(nj + 712)]. In this manner, one can use optical fields having wavelength A to 
produce atom gratings having periodicity much less than A. 
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I. INTRODUCTION 



Over the past several years, atom interferometry has 
emerged as an important new research area in atomic, 
molecular and optical physics The technology has 
improved to the point where it is now possible to control 
the center-of-mass motion of atoms using either micro- 
fabricated gratings || or optical fields ||. Accompany- 
ing the developments in atom interferometry and atom 
optics have been attempts to produce nanostructures us- 
ing atom optics elements. Probably the most successful 
method to date employs standing wave optical fields to 
focus atoms to a series of lines or dots having size on the 
order of tens of nanometers Q] . The period of the struc- 
tures produced in these focusing schemes is A/2, where 
A is the wavelength associated with the standing wave 
fields. Spatial features having dimensions as small as 
A/8 have been achieved by exploiting the ground state 
optical potentials for a magnetically degenerate ground 
state ||, but the period of the entire pattern remains 
equal to A/2. 

We have described previously a method for produc- 
ing spatially modulated atomic densities 0. When an 
atomic beam passes through one or more standing wave 
optical fields having wavelength A, it is possible to use co- 
herent transient techniques to create atomic "gratings" 
having spatial period equal to A/(2n), where n is a posi- 
tive integer. The atomic gratings arise as a result of the 
nonlinear interaction between the atoms and the fields. 
For example, when an atomic beam passes through a res- 
onant, standing wave optical field, the field can create all 
even spatial harmonics in the excited and ground state 
populations. As a result of spontaneous emission, the 
excited state gratings decay back to the ground state; 
however, for properly chosen level schemes the ground 
state gratings are not completely "filled in." As a con- 
sequence, one has long-lived ground state gratings with 
which to operate. The fate of these ground state gratings 
depends on the collimation of the atomic beam. 

If the angular divergence of the atomic beam is 9b, the 
gratings persist for a distance of order Lb = A/(2n0&) 



following the atom-field interaction. For distances larger 
than Lb, the gratings wash out quickly as a result of 
Doppler dephasing. The grating structure is not lost, 
however. If the atoms interact with a second standing 
wave field, the various spatial harmonics will rephase at 
different distances following the interaction with the sec- 
ond field. By a proper choice of optical field strengths, 
one can create and isolate atomic density patterns that 
correspond very closely to pure, higher order spatial har- 
monics pj. The patterns could be deposited on a sub- 
strate and used as elements in soft x-ray systems. 

One can also use nonresonant optical fields to interact 
with the atoms. In this case, one creates a ground state 
phase grating rather than an amplitude grating j7| . If a 
highly-collimated atomic beam is used, the phase grating 
can focus the atoms to a series of lines having spatial pe- 
riod equal to A/2. For a beam having angular divergence, 
echo techniques can be used to rephase and isolate the 
various spatial harmonics [^|j8| . 

If the goal of an experiment is to create a pure, higher 
order spatial harmonic, it would be helpful to eliminate 
the lower order harmonics from the outset. In this paper, 
we describe a method in which a single field interaction 
zone can be used to produce gratings having spatial pe- 
riod equal to A/ (2n) where n > 2 (see Fig. [j]). An atomic 
beam passes through an interaction region in which three 
fields act. One of the fields has frequency fl and propa- 
gates in the z direction, while two additional fields, each 
counterpropagating relative to the first field, have fre- 
quencies Q + Si and Q + 82, respectively. A two-photon 
process in which a photon of frequency £1 is absorbed and 
one of frequency Q + 5\ or fi + 82 is emitted (see Fig. ||a) 
produces a contribution to the ground state amplitude 
varying as e 2lkz , k = 2ir/\. Such two-photon processes 
are not resonant, however, and lead to a vanishingly small 
second harmonic amplitude if \SjT\ 3> 1 (j = 1,2), where 
T is the pulse duration in the atomic rest frame. On the 
other hand, an elementary four-photon process involving 
the absorption of two photons of frequency f2 and the 
emission of one each of frequency f2 + <5i and f2 + 62 is 
resonant, provided Si = —82 (see Fig. ||b). This process 
is responsible for the creation of the fourth spatial har- 
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monic in the ground state amplitude. Thus the geometry 
of Fig. ||b can be used to produce a ground state ampli- 
tude having spatial period A/4 provided Si = —S 2 - With 
other choices of Si and 82, one can suppress an arbitrary 
number of lower order harmonics and produce a ground 
state amplitude having period A/(2n) where n > 2. 
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FIG. 1. Atom- field geometry. All fields overlap in the in- 
teraction region. The detunings Si and 82 are chosen such 
that \8\ I T, |^2 1 T 2> 1, where T is the interaction time in the 
atomic rest frame. 
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The suppression of lower order harmonics has impor- 
tant implications for nanolithography. It is shown below 
that the field amplitudes and detunings can be chosen 
such that a single, higher order atomic grating having 
period A/(2n) can be created to a high degree of accu- 
racy. Both amplitude and phase gratings can be pro- 
duced. Phase gratings can lead to atom focusing into 
lines separated by A/(2n) <C A/ 2 for n ^> 1. Amplitude 
gratings represent a pure, higher order spatial harmonic 
that can be deposited on a substrate or used directly to 
scatter soft x-rays. 

Of course, there are other methods for suppressing 
lower order spatial harmonics. A coherent beam split- 
ter can produce momenta separation of state amplitudes 
that are greater than 2hk. As long as these momenta 
components are not spatially separated, the atomic beam 
contains density matrix elements in which lower order 
spatial harmonics are suppressed. Beam splitters based 
on higher order Bragg scattering || can lead to ground 
state densities corresponding very closely to pure, higher 
order spatial harmonics. Adiabatic rapid passage ||To[| can 
also lead to very nearly pure, higher order spatial har- 
monics, albeit between state amplitudes corresponding 
to different internal states. Beam splitters based on tri- 
angular optical potentials such as magneto-optical beam 
splitters (TTJ and bichromatic beam splitters |12j produce 
large momenta separations, but not pure spatial har- 
monics. Our approach is most closely related to that of 
Yablonovich and Vrijen [|f| , who use different frequency 
fields to increase spatial modulation in two-photon mi- 
croscopy. Here we extend and expand their ideas to the 
domain of atom optics and atom interferometry. 

The paper is organized as follows: In Section II, an 
illustrative example is considered in which the second or- 
der spatial harmonic is suppressed. Methods for creating 
both amplitude and phase gratings are discussed, as are 
the effects of spontaneous decay during the pulse. Both 
numerical solutions and an approximate analytical solu- 
tion is considered in this section. In Sec. Ill, numerical 
solutions are presented for the suppression of harmon- 
ics beyond second order. It is shown that high contrast, 
amplitude gratings can be created for arbitrary n. A 
summary of the results is given in Sec. IV, and the im- 
plications for nanolithography are discussed. 




FIG. 2. Elementary processes involving fields Eo, Ei and 
E2. (a) Two photon transitions involving field _Eo and either 
of fields Ei or E2 are not resonant (the diagrams are drawn 
with Q = uj). (b) A four photon transition is resonant pro- 
vided that Si — —82- 



II. ILLUSTRATIVE EXAMPLE 



Many of the features of harmonic suppression can be 
illustrated for the field geometry of Fig. [l| The atomic 
beam propagates in the x direction and the fields prop- 
agate along the z-axis. The total field can be written as 
E(r, t) = E(z,t)f(x,y), where f(x,y) is a spatial mode 
function, 



E(M) = ~y 
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E 9 e 



i[-kz—(p+8 3 )t] 



+ C.C. 



(1) 



and c.c. stands for complex conjugate. Fields 1 and 2 
propagate in the -z direction and field in the +z direc- 
tion. In this section, we set E\ = E2, take E\ and Eq to 
be real, and set 5i = —62 = 5, such that 



-iQt 



[E e lkz 



2E x e~ ik " o.s 



(St)] + c.c. (2) 



In the atomic rest frame these fields appear as a radia- 
tion pulse and the field amplitudes become functions of 
time. The atoms are modeled as having two levels, 1 and 
2, separated in frequency by to. The problem divides into 
two parts, interaction of the atoms with the fields and 
free evolution of the atoms following the field interac- 
tion. In the main body of the paper we consider only the 
field interaction region. In the field interaction region, all 
effects associated with quantization of the atoms' center- 
of-mass motion, as well as any transverse Doppler shifts, 
are neglected. In Sec. IV, we will discuss the free evo- 
lution of the atoms following their interaction with the 
field. 

In the atomic rest frame, the state amplitudes a\ and 
d2 evolve as 



a = iF(z, t)a — 2iG(z, t)a, 



(3) 



where 



«i 

«2 



F(z,t) 
G(z,t) 



A/2 Xo (z,t)* 



cos(8t), 



(4) 



Xo (z,t) -A/2 

X i{z,t) 
X i(z,t) 

Xo{z,t) and xi(z,i) are Rabi frequencies defined by 

Xj (z,t)= Xj (t)e ik ^; 

Xj (t) = -liEj(t)/2h = [ Xj (t)]* > 0; {3 = 0,1,2}; 
ko = —ki = -k 2 = k, (5) 

A = Q — lj is an atom-field detuning, fi is a dipole mo- 
ment matrix element, and Ej(t) is a pulse envelope func- 
tion in the atomic rest frame. Spontaneous decay during 
the pulse has been neglected for the moment. Equation 
(H) can be solved numerically for any pulse envelopes 
Ej(t). The advantage of the choice of parameters con- 
sidered in this section is that they allow for a very good 
approximate analytical solution that illustrates the rele- 
vant physical concepts. 

Since the two matrix elements of G have the same time 
dependence, it is convenient to write a solution to Eq. (j^) 
in the form 



a(t) = S(z,t) a(t); 
S(z,t) = exp(-2i / G(z,t')dt' 



(6a) 
(6b) 



where the matrix a(t) satisfies the differential equation 

da/dt = -iS f (z, t)F(z, t)S(z, t)a. (7) 

For a smoothly varying pulse having duration T 3> 8 , 
the integral in Eq. (|6h| ) can be evaluated asymptotically. 
It follows that 



S(z, t) ~ exp 



-2i 



.s'm(St) 



xi (z,ty 

Xi{z,t) 

ikz 



where 



= 2 X i(t)sm(5t)/5, 



(8) 



(9) 



1 is the unit matrix, and we have used the fact that 
Xi(z,t)/ \xx(z,t)\ = exp(-ikz). The S-matrix reduces to 
the unit matrix as t ~ ±00. Combining Eqs. (^)- (|^), 
one finds 



da/dt 



A 

-i— cos(26») 



sin(26») 



1 
-1 

iV 
Ae~ ikz 



-Ae ikz 
N* 



M 
M* 

a 



where 



M = xo{t) [cos 2 (9)e- tkz + sm 2 {0)e 3lkz ] 
N = 2i Xo (t) sin(2/cz). 



a (10) 
(11) 



(12a) 
(12b) 



Equation @ can be simplified considerably if 

5>Xo(i), A, xo(t)xl(t)/S 2 - (13) 

In that limit, it is possible to "course-grain" Eq. ( ^o|) on 
a time scale greater than (5 _1 and replace all trigonomet- 
ric functions appearing in Eqs. (|l(j) and ([l2|) by their 
time averages. Using Bessel function expansions for the 
trigonometric functions, one finds 

i + Jb[4xi(t)/ai 



a = cos 2 ((9) 



p = waT(B) 



l-J [A X i(t)/8] 



sin(26») = 0, (14) 



where Jo is a Bessel function and the bar indicates a time 
average. In this limit, Eq. (fill) is replaced by 



da/dt =-i^J [i X i(t)/S] (l _° 



ae 

ae lkz + (3e- 3lkz 



-ikz 



(3e 



3 ikz 







a. (15) 



Equation ( |T5| ) admits solutions which represent both am- 
plitude and phase gratings in the ground state atomic 
density. We examine these cases separately. 
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A. Amplitude gratings 



A = 



It is possible to obtain amplitude gratings having max- 
imum contrast for the ground and excited state popula- 
tions by choosing A = (f2 = to). In this limit, Eq. (IT" 
reduces to 



da/dt = —iXo{t) 







-3ikz 



■/3e 



'Mkz 



(16) 



It will prove useful to look at the perturbative limit of 
Eq. (|l6|). When Xi(t)/8 <C 1, the time evolution of state 
2 is 



da 2 /dt=-i Xo (ty k *a 1 -2iX° it)b(l{t)Y 



S 2 



-3ikz ~ 



(17) 



The fundamental processes responsible for excitation to 
state 2 from state 1 are shown schematically in Fig. 0. 
There can be direct, resonant excitation to state 2 by 
field Eq and three-photon, resonant excitation involv- 
ing the absorption of one photon each from fields E\ 
and E2 and emission of one photon into field Eq (as 
noted above, the amplitude modulated field used in this 
section can be considered as a sum of two fields hav- 
ing frequencies f2i = Q + 81 and Q2 = ^ + 82, with 

51 = —82 = 8). The three-photon processes are resonant 
since ^2 — ^ = w when Q = u). In these diagrams, 
the Rabi frequency xo is a shorthand notation for xo( z i t) 
and contains a factor e %kz and, similarly, Rabi frequen- 
cies Xi and X2 contain factors e~ lkz . The overall ampli- 
tude for the one photon process varies as Xo^ lkz and as 
(XoXiX2) e ~ 3lkz / for the three-photon processes. The 

5 2 factor reflects the fact that the two intermediate states 
in the three-photon processes are each off-resonance by 
an energy h~S. In taking the square of the amplitude, one 
finds a spatial modulation at Akz. As seen in Sec. Ill, di- 
agrams of this type can be used to estimate the required 
field strengths for suppression of higher order harmonics. 



Xo 




FIG. 3. Elementary processes leading to the resonant ex- 
citation of state 2. The Rabi frequencies \j are a shorthand 
notation for \j{z,i) and contain spatial phase factors. In 
forming the upper state probability, interference terms lead 
to a spatial modulation having period A/4. 



Equation (|16|) cannot be solved analytically for arbi- 
trary pulse envelope functions. On the other hand, it 
can be solved analytically if we assume square pulses, 
Ej(t) = Ej for < t < T and zero otherwise. [Equation 
(0) remains valid for a square pulse. Moreover, if one 



chooses a detuning ST — nn, then the S-matrix in Eq. 
@ reduces to the unit matrix at the beginning and end 
of the pulse. Thus one can use Eq. (|l6|), even though a 
square pulse does not satisfy the adiabaticity conditions 
at turn on and turn off.] For a square pulse, one finds 

MT)| 2 = \~a 2 (T)\ 2 = sin 2 [ Xo T \ae^ kz + (3e 3lkz |] 
= sin 2 {xoT([(l + </ 2 {4xi/£}) 

+ (1 - J 2 {4xi/<5}) cos4fcz] /2) 1/2 } . (18) 

It follows from Eq. ( |l8| ) [or from Eq. ( |T6| )] that the ex- 
cited state population immediately after the pulse is a pe- 
riodic function of z having period A/4. The second order 
spatial harmonic having period A/2 has been suppressed, 
owing to the large detunings ±<5 of the fields propagat- 
ing in the — z direction. By taking a field strength corre- 
sponding to the first zero of the Bessel function Jo [4xi /8] 
(Xi ~ 0.65), one can optimize the grating contrast for the 
smallest possible value of XoT. If Jo[4xi/6] = 0, then 
|a2(r)| 2 = sin 2 [xoT cos(2kz)]. This function is plotted 
as a function of kz in Fig. || for several values of XaT- 
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FIG. 4. Graphs of the upper state probability as a function 
of kz for excitation by a rectangular pulse, with Si — —82 = 8. 

For values of XoT <C 7r/2, the grating contrast in the 
excited state amplitude is less than unity, but the grat- 
ing is close to pure fourth harmonic. For xqT > 7r/2, the 
grating contrast is always unity, but harmonics higher 
than fourth are evident. A choice of pulse area that 
leads to high grating contrast while producing a spa- 
tial distribution very close to pure fourth harmonic is 
XoT = 1.3, for which the Bessel function expansion of 
sin 2 [xoT cos(2kz)] yields 

|a 2 (T)| 2 - [1 - J (2.6) + J 2 (2.6)cos(4fcz) 
- J 4 (2.6) cos(8fo) + . . .] /2 
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w [1.1 + 0.92 cos(4/cz) - 0.17cos(8/cz)]/2. (19) 

This grating has a fringe contrast of 0.93 and consists 
primarily of fourth harmonic. 

The qualitative results do not change when smooth 
pulses are used. The original differential equation (^) 
was solved numerically for Gaussian pulses, 

X J (<) = X,e- (i/T)2 /^- (20) 

From this point onward, we set T = 1 and evaluate all 
times in units of T and all frequencies in units of T _1 . 
In these units, the quantity 2\j corresponds to the area 
of pulse j. In all the examples, 5 > 3, which, for the 
Gaussian pulse (|20|) , is sufficient to guarantee adiabatic- 
ity, provided that conditions ( |l3| ) are also met. 

In Fig. ||, the solutions of the exact equations [Eq. (^)] 
for the upper state probability are compared with those 
of the coarse-grained equations [Eq. (jlq)]. Recall that 
the solutions of the course grained equations should ap- 
proach those of the exact equations as the ratio S/xo in- 
creases. This feature is seen clearly in Fig. 5. The curve 
having S/xo = 40 correspond to the solution of both the 
exact and course-grained equations - they are not distin- 
guishable for this ratio of S/xo- Even for 5/xo = 5> the 
results do not differ by much. 



and xi> moreover, the solution has period A/2 rather 
than A/4. Although not evident from the figure, the 
<5 = Xo results differ from pure A/4 periodicity by as 
much as 10%. These differences are more evident in Fig. 
||, drawn for 5/xo = 1/3- To assure A/4 periodicity, one 
must have S/xo ^ 5. 




0.(10 



Jl/2 It 

kz 



Gaussian Pulse x l ,=" X|=x,=(l2-5/l5)8 A=0 




FIG. 5. Graphs of the upper state probability as a function 
of kz for excitation by a Gaussian pulse, with 8i = —82 = 8. 
In this and subsequent graphs, all frequencies are in units of 
T -1 . With increasing ratio of S/xo, the solution of the exact 
equations (Q) approaches that of the course-grained equations 
(p_7[), represented by the solid line in the plot. 

The course grained equations depend only on the pa- 
rameters xo and <Vxi; consequently the solutions of the 
exact equations depend on these two parameters only for 
d/Xo ^ !• If <V%o ^ 1) the solution of the exact equa- 
tions depends independently on the parameters xo, 



FIG. 6. Graphs of the upper state probability as a function 
of kz for excitation by a Gaussian pulse, with Si = —62 = S 
and 5/xo = 1/3. For this ratio, the "course-graining" approx- 
imation is no longer valid and the period of the amplitude 
grating is A/2 rather than A/4. 
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FIG. 7. Graphs of the upper state probability as a function 
of kz for excitation by a Gaussian pulse, with 81 = —62 = S 
and several values of xo- 
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In Fig. |t], the upper state probability is plotted as a 
function of kz for fixed Xi/<5 = 5/6, and several values of 
Xo- These curves follow the same qualitative behavior as 
those shown in Fig. || for the square pulse. In Fig. [| the 
upper state probability is plotted as a function of kz for 
fixed xo = t and several values of Xi/<5- For xi/£ < 0.5, 
the grating contrast is less than unity, but the grating is 
very nearly pure fourth harmonic. With increasing Xl/<^ 
the grating contrast approaches unity and harmonics of 
order (4 + 4n), n > 1, are evident. 



plitude, it is convenient to use the semiclassical dressed 
states associated with Eq. (|l6|). These states are ob- 
tained by instantaneously diagonalizing Eq. ([l6|). If 



a2Jq2 /4xiW 



+2xo(*) 2 



1- J 2 



(l + Jo < 


IT}) 


/4xiWl 


>^ cos4fczj 


I « J 



> 1 



(21) 
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FIG. 8. Graphs of the upper state probability as a function 
of kz for excitation by a Gaussian pulse, with Si = —82 = 5 
and several values of xi- 



the system remains in the instantaneous eigenstate that 
evolves from the ground state as the field is turned on 
and returns to the ground state following the pulse. The 
net phase change is simply the integral of the dressed 
state energy divided by Ti. Explicitly, one finds. 



ai(oo) 



(22) 



where 



dt { A 2 Jg 



+2xoW' 



1- J A 



1 + J, 



4xi (Q 
S 

4X1 (*) 



4X1 (*) 



6 

cos Akz 



1/2 



(23) 



For Gaussian puls es (|20| ) , the phase </> is plotted as a func- 
tion of fcz in Fig. |10| for several values of xo, Xi> an d 
A. 



B. Phase Gratings A > 1 

When A > 1, the elementary processes shown in Fig. || 
are no longer resonant. Thus one expects that the excited 
state population following the pulse to be negligible. On 
the other hand the resonant, elementary 4-photon pro- 
cesses, shown in Fig. ||, lead to a modification of the 
phase of the ground state amplitude. 




FIG. 9. Elementary processes leading to a phase grating in 
the ground state amplitude. The second and third diagrams 
in the figure are spatially modulated with period A/4. 

To calculate the phase change of the ground state am- 
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FIG. 10. Graphs of the ground state phase as a function 
of kz for excitation by a Gaussian pulse, with Si = — 82 = 8. 
The solid lines are a solution of the exact equations while the 
dashed lines are obtained from the course-grained, dressed 
state equations. 
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The solid line represents solutions of the exact equa- 
tions (0) and the dotted line the values given by Eq. 
(p3|), with the phase arbitrarily taken equal to zero at 
kz = 0. The results agree when conditions ( |2l| ) and (|i~3| ) 
are satisfied. For %i/<5 > 0.6, there are times during the 

atom-field interaction for which Jo j 4x * ^ | = 0. As a 

consequence, condition ( pl| ) is violated near kz = it/ 4 
and 37r/4. This explains the small deviations between 
the exact and approximate solutions for these values of 
kz. For all cases shown, |ai(oo)| ~ 1. It is seen that 
significant phase modulation of the ground state can be 
achieved with suppression of the second spatial harmonic. 

A necessary condition for significant spatial modula- 
tion of the phase is Xoxl/6 2 > 1. If A 2 J 2 {%^} < 

Xo(t) 2 , then the effective pulse area for the creation of 
the phase grating is of order Xoi while it is of order Xg/A 

for A 2 J 2 {%W}» Xo (t) 2 . 



C. Role of Spontaneous Emission 

So far we have neglected the role of spontaneous emis- 
sion during the excitation pulse. To investigate the role of 
spontaneous emission during the pulse, we adopt a highly 
simplified decay scheme in which state 2 decays to a state 
outside the two-state manifold with rate 27. As a result, 
population leaks out of the two-level manifold. Clearly, 
any excited state amplitude grating will be diminished 
by decay. It also turns out that phase gratings are de- 
stroyed by spontaneous decay since the excited state pop- 
ulation during the pulse is not negligible under conditions 
favoring a significant spatially modulated phase grating. 
Spontaneous emission breaks the adiabatic following that 
would return the atom to its ground state in the absence 
of such decay . 

Within the context of this model, spontaneous emis- 
sion can be included in Eq. (0) by the addition of a term 




1 



This decay term results in an addition to Eq. (|T^) of the 
form 



l-Jo[4 X i(t)/<5] 
2 







1+Jq[4xi(*)/5] 
2 



For any nonvanishing field strength xii both dressed state 
amplitudes a\ and 5,2 decay. For the field strength giving 
rise to maximum spatial modulation, Jo[Axi(t)/S] ~ 0, 
spontaneous emission destroys grating formation when 
7 > 1. This is easily understood. To have strong spa- 
tial modulation of amplitude or phase gratings, there is 
always a significant excited state population during the 
pulse; any decay results in a loss of population from the 
system in our simple decay scheme. 



We have also carried out density matrix calculations 
using various decay schemes. It appears that one must 
require 7 < 0.5 to retain reasonable grating contrast. 



III. SUPPRESSION OF HIGHER ORDER 
HARMONICS 

One can generalize the previous results to allow for the 
suppression of harmonics higher than second order. The 
appropriate field geometry is shown in Fig. Sponta- 
neous emission during the atom-pulse interaction is ne- 
glected. For the field ([!]), the state amplitudes evolve 
as 



a = — iF(z, t)a — iH(z, t)a 
where F(z,t) is defined by ([J) and 



(24) 



H(*,t) 



( 

X i(z,t)e- iS ^ 
V +X2(z,t)e~^ t 



(25) 



Since H(z, t) cannot be factored as a constant matrix 
times a function of t, the method used in Sec. II is not 
applicable. On the other hand it is not difficult to obtain 
numerical solutions of Eq. ( ^5|) . 

Diagrams similar to those shown in Figs. || and ^| can 
help us to estimate the field strengths needed to obtain 
both amplitude and phase gratings. Before considering 
such diagrams, we note that the lowest order ac Stark 
shift, 



Z={xl/h + xllh) 



(26) 



no longer vanishes automatically as it does for an ampli- 
tude modulated field {xi = X2', S\ = —62)- To suppress 
the higher order harmonics, it is necessary that |<5i| and 
I #2 1 be large compared with all characteristic frequencies 
in the problems, including the ac Stark shifts. In order 
to have |<Si|,|<$2| > h is necessary that xi/^i ^ 1> 
X2/ ^2 ^ 1j unless 5i and S2 have opposite signs and the 
field strengths are chosen such that 



(xi/Si + xt/52) = 0. 



(27) 



In this case it is possible to have xi ~ ^1 an( i still have 
l^i I j 1^1 > £■ Since it is necessary to have xi ~ ^1 f° r 
significant spatial modulation in the state populations, 
it is important to choose the Rabi frequencies such that 
Eq. (|27| ) is satisfied. Higher order contributions to the 
ac Stark effect may still play a role. From the numer- 
ical solutions, it appears that one must choose xi anc l 
X2 less than or of order |<5i|,|#2| > to obtain the desired 
periodicity. 

The elementary processes giving rise to an amplitude 
grating are illustrated in Fig. [ll]. The grating is formed 
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by terms representing the interference of the two dia- 
grams shown in the figure. 



a 9 ~ xo 



i- 



+ Xi 




FIG. 11. Elementary processes leading to the resonant 
excitation of state 2 and suppression of spatial harmonics 
higher than fourth. The detunings are chosen such that 
niSi — 0, leading to an amplitude grating having pe- 

riod A/[2(m +n 2 )]. 



Condition @ can be satisfied if xo ~ p-Vim+m), 
but, to suppress lower harmonics, it is also necessary that 
Xo <C S. As (ni + n 2 ) increases, S and xo must be taken 
larger and larger to simultaneously satisfy these condi- 
tions. 

In Fig. |l2| the upper state population is plotted as a 
function of kz for {n\ = 2, n% = 1) (sixth harmonic) 
and for (n\ = 3, n 2 = 2) (tenth harmonic). The field 
strengths are chosen such that \2 — \J n i/ n 2Xi accord- 
ing to Eqs. (27) and (p8|), and \i is taken to be of order 
5\. It is seen that amplitude gratings having periods as 
small as A/10 can be produced with contrast of order 
unity. For sufficiently large 5, an arbitrary number of 
lower order harmonics can be suppressed. 



The resonance condition for the second diagram is 



ni<5i + n 2 5 2 = 



(28) 



when A = and the integers n% > n 2 contain no common 
factors. The lowest order nonvanishing harmonic varies 
as cos[2(ni + n 2 )kz}. The amplitude associated with this 
harmonic is of order 



A(m,n 2 ) 



Hi -\-no n-i no 

Xo Xi X 2 
p(m+n 2 -i) 



(29) 



where S is of order |<5i| or \S 2 \. To obtain good grating 
contrast, one requires that 



+ Xi 



1- 




FIG. 13. Elementary processes leading to a ground 
state phase grating and suppression of spatial harmonics 
higher than fourth. The detunings are chosen such that. 
niSi + U2&2 = 0, leading to a phase grating having period 
A/[2(m +n 2 )]. 



A(m, n 2 ) ~ 1. 



(30) 



Gaussian Pulse A=0 



X„=5t, X,=29, x,=2W 2, S,=30, 8 2 = -60 

!»' r Xq=9», X|=50, x]=50V(3ffl), 5,=50, -75 




FIG. 12. Graphs of the upper state probability as a 
function of kz for excitation by a Gaussian pulse, with 
ni8\ + U2&2 = 0. Graphs are shown for (m = 2, ri2 = 1) 
(sixth harmonic) and (m = 3, ri2 = 2) (tenth harmonic). 



Gaussian Pulse A=5 



■ Xo= 20, X,=50, X 2 =5(W 2, 8,=50, 8,= -100 
X =125. x,=40o"x,=400V(3/2), 8,=400, S 2 = -600 




FIG. 14. Graphs of the ground state phase as a function of 
kz for excitation by a Gaussian pulse, with ni5\ + U2&2 = 0. 
Graphs are shown for (ni = 2, 712 = 1) (sixth harmonic) and 
(ni = 3, ri2 — 2) (tenth harmonic). 

The conditions for producing phase gratings are simi- 
lar, but a little more severe, than those for producing am- 
plitude gratings, provided one takes A > 1, but xo A. 
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The elementary process es g iving rise to an phase grating 
are illustrated in Fig. [13|. The highest contrast phase 
gratings are produced when A(nx,ri2) ~ A, which can 
be satisfied if xo ~ ( y 1 - 2 /(«i+«=)A 1 /( ni+ " 2 ). Along with 
the requirement that xo 8, this condition implies that 
large values of S and xo are needed to produce large 
phase gratings as (ni + n^) increases. The phase of the 
ground state amplitude is shown as a function of kz in 
Fig. |l4| for (n\ — 2, m — 1) (sixth harmonic) and for 
(ni = 3, ri2 = 2) (tenth harmonic). In all cases shown, 
|ai(oo)| ~ 1 . Clearly it is also possible to produce phase 
gratings having reasonable contrast with the suppression 
of lower order harmonics. 



IV. DISCUSSION 

It has been shown that atom-field interactions can 
be used to create high contrast amplitude gratings for 
atomic state populations and phase gratings for atomic 
ground state amplitudes. If countcrpropagating fields 
having wavelength A are used, it is possible to choose 
atom-field detunings such that the lowest spatial har- 
monic component in the gratings has period A/2n, where 
n is an integer. For the sake of definiteness, we will refer 
to the harmonic having period A/2n as the (2n)th har- 
monic. The gratings produced are not pure in the sense 
that they contain all integer multiples of the (2n)th har- 
monic, but the field strengths can be chosen in such a 
manner as to maximize the contribution from the (2n)th 
harmonic. Once the gratings are created, the question 
remains as how to image the gratings at some distance 
L from the atom-field interaction zone. One can probe 
the gratings by applying counterpropagating fields hav- 
ing frequencies f2 and Q+Si to generate a signal with fre- 
quency Q + 62 ; alternatively, one can deposit the atomic 
density grating on a substrate. The situation is different 
for amplitude and phase gratings. 

In the case of amplitude gratings, one has the op- 
tion of working with distances L that correspond to ei- 
ther classical or quantum-mechanical scattering for the 
atoms. If L > Lt, where the so-called Talbot length Lt 
is defined by Lt = 2{X/2n) 2 / XdB (A^s is the atomic de 
Broglie wavelength), the atomic center-of-mass motion 
for the atoms must be treated quantum mechanically; 
for L <C Lt, this motion can be treated classically. Typ- 
ically, Lt is of order of a few cm for a thermal beam. We 
restrict our discussion to the classical scattering limit. 
For classical scattering, the total atomic density is not 
modulated, as long as no state-selective mechanism re- 
moves atoms from the system. In other words, the ex- 
cited and ground state population can be modulated fol- 
lowing interaction with the fields, but the sum of these 
populations is not modified. If the excited state popu- 
lation decays to the ground state following interaction 
with the field, the grating structure disappears. There 
are several methods to maintain the gratings following 



the atom-field interaction. The most direct method is to 
use optical transitions, such as those in Ca, Mg, or Yb, 
with lifetimes approaching a microsecond or longer. One 
can then probe the population difference grating follow- 
ing the atom-field interaction or selectively ionize excited 
state atoms to leave a net atomic density grating that 
could be deposited on a substrate. 

For shorter-lived atomic transitions, there are several 
alternatives. One possibility is to ionize some of the ex- 
cited state atoms during the atom-field interaction, which 
would leave a net atomic density grating. Another pos- 
sibility is to use atoms whose ground state consists of a 
manifold of substates. When the excited state decays, 
the total ground state density is not modulated, but that 
of specific substates may be modulated. By selectively 
ionizing some of the ground state sublevels, one again 
achieves a net atomic density grating. Still another al- 
ternative is to use a beam of metastable atoms and drive 
transitions between the metastable state and an excited 
state that can decay back to the ground state p5| . Fol- 
lowing the atom-field interaction, there will be a grating 
of metastable atoms that can be detected selectively. 

In the case of shorter-lived excited states, one can run 
into problems related to the effective pulse interaction 
time T which corresponds to the transit time of the atoms 
through the field interaction zone. Recall that all fre- 
quencies are measured in units of T _1 . For appreciable 
gratings, it is necessary that 7T < 0.5. If the fields are 
focused to a spot size of 10 \x and the longitudinal atomic 
speed is 5 x 10 4 cm/s, then T w 20 ns. For these param- 
eters, one must use atoms whose excited state lifetime 
exceeds 20 ns. It is also possible to limit the atom-field 
interaction time by using pulsed fields. For transit times 
of order 20 ns, the detunings and Rabi field strengths 
considered in this paper correspond to frequencies in the 
MHz to GHz range. Such values do not pose any serious 
experimental difficulties. 

Assuming that we can produce a density grating, many 
methods are available for imaging the gratings. Perhaps 
the most direct method is to detect the grating immedi- 
ately following the atom-field interaction. This scheme is 
analogous to free-induction decay (FID) in coherent tran- 
sient spectroscopy. All harmonics in the beam are super- 
imposed, and the grating begins to wash out as a result of 
the transverse Doppler effect at distances L(, = A/ (2n6f,), 
where 6^ is the atomic beam divergence. To avoid overlap 
of the different harmonic components, echo techniques 
can be used With two interaction zones separated 
by a distance Lq, different harmonics are focused at dif- 
ferent distances following the second interaction zone. 
For example, one can optimize the (2n)th harmonic at 
L = 2Lq, the (4n)th harmonic at L = ZLq/2, the (6n)th 
harmonic at L = ALq/3, etc. In this manner one is able 
to isolate and focus harmonics higher than 2n which were 
created in the first field interaction zone q|. For exam- 
ple, one could create 4th harmonic in the first interaction 
zone and 16th harmonic in the second interaction zone 
to optimize the 12th harmonic at L = ALq/3. 
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Although FID and echo techniques can be used also 
for phase gratings § , we consider here only the focusing 
capabilities of phase gratings. The phase grating can be 
thought to arise as the result of an ac Stark shift poten- 
tial that acts as a spatially modulated index of refraction 
for the atoms and causes them to focus at a distance of 
order Lt/A, where A is some effective pulse area in the 
problem associated with the creation of the (2n)th har- 
monic. In this manner, one can focus an atomic beam 
to a series of lines (or dots if a 2-D geometry is utilized) 
having spacing X/2n. 

In summary, we have outlined methods for generating 
and detecting spatially modulated atomic density ampli- 
tude or phase gratings having period X/2n using optical 
fields having wavelength A. All harmonics having larger 
periods are suppressed by a proper choice of atom-field 
detunings. 
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